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Abstract^ 

The problem of the determination of the Hilbert-space metric O which ren- 
^ I ders a given Hamiltonian H self-adjoint is addressed from the point of view 

of applicability of computer-assisted algebraic manipulations. An exactly 
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solvable example of the so called Gegenbauerian quantum-lattice oscillator 
is recalled for the purpose. Both the construction of suitable = Q{H) 
(basically, the solution of the Dieudonne's operator equation) and the deter- 
O I mination of its domain of positivity are shown facilitated by the symbolic 

algebraic manipulations and by MAPLE-supported numerics and graphics. 
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(N 1 Introduction 

in 



In the series of papers [I] - [7] we interpreted the Dieudonne's quasi-Hermiticity 
relation [S] 

H^Q = QH (1) 



as an operator equation which connects a given, "input" quantum Hamil- 
tonian H with an unknown, "output" operator called the Hilbert-space 
metric of the quantum system in question ([9j; cf. also a more detailed ex- 
planation of the underlying physics in Appendix A below). 

In all of our papers we felt addressed by the underlying physics (i.e., by 
quantum mechanics in its form described, say, in Refs. [TOl E]) and did not 
pay too much attention to the underlying constructive mathematics. In our 
present paper we intend to fill this gap by redirecting our attention to the 
computer-assisted symbolic-manipulation background of our results. 

We shall assume that both of the operators H and B in Eq. ([1]) are defined 
in a vector space V with the Dirac-ket elements {ip) G V. For the sake of 
simplicity we shall further assume that dimV = < oo. In such a setting 
one may select various toy-model matrices H^^^ and construct the eligible 
metrics 0. For such a purpose we shall choose the Gegenbauerian quantum 



^Work supported by the GACR grant Nr. P203/11/1433, by the MSMT "Doppler 
Institute" project Nr. LC06002 and by the Inst. Res. Plan AV0Z10480505. To be 
presented to the Computer Algebra in Scientific Computing (CASC 2011) conference in 
Kassel, Germany, during September 5-9, 2011. 



A^-site lattices and study H = H'^^^a) and = Q^^\a) of Ref. [2] (cf. 
section 12]). The reasons of a facilitated algebraic tractability of Eq. ([1]) will 
be then clarified in section [3l with some complementary numerical aspects 
mentioned in section |H 



2 Gegenbauerian quantum A^— site lattices 



The general Hermitian conjugation prescribed by Eq. (ITSl) of Appendix A 
must be compatible with the principles of Quantum Mechanics. This means 
that, formally speaking [TT], our choice of the metric G must guarantee the 
Hermiticity of the observables (i.e., in our present paper, just of the Hamil- 
tonian H) with respect to this conjugation [9]. Such a requirement implies 
the necessity of the validity of the above-mentioned relation ([1]). 

The latter relation will be called here, for the sake of brevity, Dieudonne 
equation. As long as this is the matrix equation, it seems to be an overdeter- 
minate constraint. Its A^^ items have to be satisfied by the mere A^(A^ + l)/2 
independent matrix elements of the general N hj N real and symmetric 
matrix B with positive eigenvalues. Via a deeper study of an illustrative 
example taken from Ref. [2] we intend to demonstrate that the situation is 
much more user friendly. 

First of all, the Hermitian conjugation may be perceived as a symmetry 
of Eq. ([1]) so that just its upper-triangular nontrivial subset remains relevant. 
This implies that the whole set of equations is in fact underdeterminate. A 
priori, the independent solutions G will form an A^— parametric family. This 
observation is compatible with the explicit constructive results published in 
Ref. [g. 

Secondly, the message delivered by Ref. [2] was aimed at the physics 
audience. Our present study will complement these results by their more 
systematic derivation and by the more explicit explanation of their formal 
structure. Keeping this aim in mind we shall consider the A^— dimensional 
matrix Schrodinger equation 



(2) 



with the prescribed bound-state eigenvectors 

/ 



(0|^f^) = G(0,a,E„) \ 
(l|V.f)) = G(l,a,E„) 



V {N-l\^iP) = G{N-l,a,E^) ) 



(3) 



where, in the notation of MAPLE [12j, the symbol G{n,a,x) denotes the 
n— th Gegenbauer polynomial G{n,a,x) equal to polynomial C'^{x) in the 
notation of Ref. [13] or to Cn\x) according to Ref. [H]. Under such an 
assumption, naturally, the explicit form of the related Hamiltonian is the 
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tridiagonal array 
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which defines the manifestly asymmetric by matrix H^^\ 

The vahdity of such an assignment is equivalent to the standard three- 
term recurrences for the Gegenbauer polynomials while the N hj N matrix 
truncation is equivalent to the implicit-equation identification of the real and 
non-degenerate spectrum a{H^^^) = {En} of the bound-state energies with 
the roots of the A^— th Gegenbauer polynomial, 

G(iV,a,E„) =0. (5) 

Naturally, such a secular equation may be considered solvable with an ar- 
bitrary numerical precision. The only nontrivial task represented by the 
complete description of the model will lie in the choice of a metric G com- 
patible with Dieudonne Eq. ([1]). The method has not been described in 
Ref. [2] because it just consisted in the brute-force insertion of a general real 
and symmetric ansatz for G(^) and in the subsequent trial and error analysis 
of Eq. (Il]) after its insertion. 

In the metric-construction problem the key difficulties are twofold. Firstly, 
the ansatz for contains too many (i.e., A^(A^ + l)/2) unknowns and there 
are no criteria for the clarification which ones of them should be selected as 
the "optimal" independent set. Even at the very small integers N, prelimi- 
nary MAPLE-based brute-force algebraic symbolic-manipulation-solution ex- 
periments starting from Eq. ([1]) and from a few randomly selected A^— plets 
of the tentative independent matrix elements of generated just the obscure 
many-page results for all of the A^— and a— dependent matrix elements of 
Q{H). The second difficulty emerged with the necessity of the parameter- 
range-specifying guarantee of the obligatorily positive-definite nature of any 
resulting A^ by A^ matrix candidate Qa{H) for the metric (characterized or 
distinguished, in general, by a suitable multiindex a). 

The core of the success (i.e., of the resolution to both of these parallel 
algebro-numerical difficulties) has been revealed to lie in an interactive and 
iterative approach to both of the problems. In more detail this approach is 
to be described in what follows. 



3 The Dieudonne's equation 

In the light of the old Dieudonne's idea [8] it seems interesting to replace 
the current textbook Hermiticity property H = of the current selfadjoint 
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Hamiltonians in quantum mechanics by the weaker assumption containing 
a nontrivial "metric" G 7^ / [9]. In this context relation ([1]) guarantees the 
reahty of the energies provided only that we require that the operator 9 = 6^ 
is, roughly speaking [15], positive and invertible, i.e., tractable as a metric 
in the Hilbert space of states where the superscript means "standard" . 

The recent growth of popularity and applicability of the quantum models 
requiring the unusual metrics 9 = 9'^"^^ > / may be found reviewed, e.g., 
in Refs. [lOl |16l |T7] . In our present paper we circumvent a number of tech- 
nicalities by studying just the quantum models defined in finite-dimensional 
Hilbert spaces. Thus, we may identify l-L^^^ = C^. Moreover, for the sake of 
definiteness, we shall only pay attention to the models where the Hamiltonian 
matrices possess the tridiagonal real-matrix form. 

3.1 An interactive algebraic-solution technique 

As an illustrative example we shall use the Gegenbauerian quantum lattice 
model of Ref. |2] described in the preceding section. In fact, in Ref. [2] we 
described the results showing the feasibility of the brute-force construction of 
the complete family of the metrics Q{H) admitted by the Dieudonne's linear 
algebraic constraints ([T]). In our present continuation of this effort we intend 
to provide a deeper insight in the problem explaining the reasons why our 
construction of the metrics appeared to be so successful. 

We have to admit that with our very specific choice of the Gegenbauerian 
model in Ref. [2] we were unexpectedly fortunate. This fact may be demon- 
strated, say, by the recollection of the similar constructive attempts based on 
a different choice of the N hj N "input" Hamiltonian as reported in Ref. 
After the construction of 9 (if) at the dimension as low as iV = 4 it has 
been argued there that the construction at the very next = 5 appeared 
almost prohibitively complicated. This is really in contrast with the results 
of Ref. |2] which proved valid at any integer N. 

The core of the dimension-independent universality of the above-mentioned 
Gegenbauerian result may be seen in the combination of the extremely sim- 
ple bidiagonal form of the Hamiltonians H^^^ (a) with the comparably simple 
a— dependence of its matrix elements. The relevance of both of these ingre- 
dients becomes obvious when we recall the explicit A^ = 4 sample of the 
Hamiltonian 
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together with the general ansatz for the metric 
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In such a setting we may study the 16-plet of the resulting relations, out 
of which Nr. 1, Nr. 6, Nr. 11 and Nr. 16 (i.e., diagonal items) remain 
trivial while the off-diagonal items form an antisymmetric matrix. Out of 
the remaining six independent items (say. Nr. 2, 3, 4, 7, 8, and 12) there 
is just one (viz.. Nr. 4) which involves just two unknown quantities (viz., 
h and c). This leads to the decision of taking d, c, b and k as independent 
parameters and of eliminating, in the first step, h via item Nr. 4, 

, c(a + l) 



2 (a + 2) a ■ 

The inspection of the new set of items reveals that the simplest one is now 
just Nr. 3 which defines g as a. function of b and d, with the next-step Nr. 8 
defining n as a function of d and (newly known) g. We are left with the three 
items Nr. 2, 7 and 12 which couple {k,f), {f,m) and {m,j), respectively. 
As long as we decided to use k as the fourth unconstrained parameter this 
means that in the same order we now define, step by step, the missing items 
/, m and, ultimately, j. The result is complete yielding 



h = 1/2 



c(a + 1) 



9 = 1/2 



a + 2) a 
ba + 3b + 2da^ + Ada + 2d 



n = 1/2 

/ = l/2 
m = 1/2 



(a + 3) a 
-Q da — 10 d + ba + 3 b 



[a + 3) a (2 a + 1) 

{2 00^ + ka + ca + 2 k){a + 1) 
(a + 2) a2 

2 ca^ + ca^ — 7ca + ka'^ + 5ka + 6k 



J = -1/4 



(a + 3) a2 (2a + l) 

6 ca^ + 10 ca — ka^ — 5ka — 6k 
a2 (2a + l) (a + 2) (a + 1) 



3.2 The case of general N 

After the above-explained heuristic exercise we are prepared to consider the 
real and symmetric general ansatz for the metric 



and prove the general result. 
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(6) 



Definition 1. At any N > 2 the insertion of the by iV Hamiltonian 
H = if (a) given by Eq. (jlj) and of the general real and symmetric matrix 
ansatz ([6]) for the metric O = Q^'^\a) defines the A^ by A^ matrix array 
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of the linear Dieudonne equations = 0. Its ordered version has the 
form Ta = with 

n = A^i.AT , 

r2 = Mi,N~i,r3 = M2,N , 
r4 = Mi,N~2 , = M2,N-i , re = M3,n , 

. . . , 

'^(Af-l)(Ar-2))/2+l = -Ml,2 , • ■ ■ , rN{N-l)/2 = Mn-1,N ■ (7) 

Theorem 1. In terms of the freely variable N—plet of the real initial param- 
eters Qij, j = 1,2, ... ,N the Dieudonne equation in its ordered version ([^ 
defines, step by step, the respective "missing" matrix elements 

02,7V-1 , 03,Af ) 
. . . , 

©2,2 , ©3,3 , ©4,4 , • • • , 0Af,iV (8) 

in recurrent manner. 

Proof. Once we revealed the diagonal-wise-arranged recurrent pattern it is 
easy and entirely straightforward to verify its validity by the corresponding 
trivial rearrangement of the two matrix multiplications in Eq. ([T]). □ 

Remark 1 . The diagonal-wise recurrent nature of Eq. ([T]) given by Theorem 
1 has only been revealed by the post factum inspection of the results of 
Ref. [2]. 

4 The positive definiteness of the metric 

During the recent years we are witnessing the remarkable growth of pop- 
ularity of the building of quantum models which combine the "false" non- 
Hermiticity H ^ W of the comparatively elementary Hamiltonian acting in 
a "friendly" Hilbert space with the simultaneous "sophisticated" Her- 
miticity H = of the same Hamiltonian in another, less usual, amended, 
"standard" Hilbert space H^^^ endowed with a nontrivial metric Q = Q^^^ ^ 
I. The essence of such an innovation becomes entirely transparent and ob- 
vious when one eliminates the partial confusion caused by the traditional 
terminology. The key point is that one never leaves the abstract theoretical 
framework of quantum theory. Just a few new mathematical tricks (like, typ- 
ically, an unusual, non-unitary generalization of the most common Fourier 
transformation) are being added to the traditional textbook recipes. 

4.1 The formulation of quantum theory using an ad 
hoc triplet of complementary Hilbert spaces 

The Dieudonne-equation constraint imposed on a Hamiltonian H is in fact 
equivalent to the manifest Hermiticity of its isospectral image 

i) = QHQ~^ = l)^ (9) 
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[TT] . In principle (though not always in practice) the latter operator is defined 
as acting in the physical Hilbert space H^^^ in which the traditional, trivial 
metric 9*^^^ = / is being used. 

Both the Hilbert spaces H^^^ and T-L^^^ may be perceived as unitary equiv- 
alent. We may deduce 

i)^ = {n-^y (10) 

After we abbreviate Q^Q := we end up with the Dieudonne's relation ([1]). 



4.2 The Gegenbauerian illustrative example 



For our real and finite-dimensional Gegenbauerian Hamiltonians H = H^^^ (a) 
which are given in advance, the Dieudonne's relation ([T]) forms the set of iV^ 
constraints imposed upon the [A^(A^-|- 1)/2]— plet of the unknown real matrix 
elements of the metric matrix = 6^ > 0. In papers |1] - |7] we proposed 
the non- numerical, symbolic-manipulation approach to constructions of a 
complete solution of this linear algebraic system. What remains for us to 
construct is the appropriate domain V of free parameters for which these 
candidates for the metric remain positive definite, i.e., truly eligible in the 
appropriate definitions of the generalized Hermitian conjugation and/or of 
the appropriate Hilbert-space inner product. 

In Ref. [2] we discussed a few specific examples of candidates Q^^\a) for 
the Gegenbauerian metrics. We revealed that such a study leads to a purely 
numerical description supporting the hypothesis that the domains T)^^\a) 
change "smoothly" with and lead to the non-empty and sufficiently large 
limiting domains limN^ooT^^^Kcb) = 'D^°°\a) ^ 0. 

A comparatively weak A^— dependence characterizes even the domains 
T)^^\a) at small N ^5. The direct evaluation of the eigenvalues of Q^^\a) 
(i.e., the more precise determination of the boundaries dV^^^a)) only suffers 
of the errors caused by the multiple-scale nature of these eigenvalues. 

In Ref. [2] we were only able to provide a transparent graphical illustration 
of the free-parameter-dependence of the spectrum of selected Q^^\a)s at 
the very first nontrivial dimension A^ = 3. In the context of programming 
in MAPLE (offering an adaptable floating-point precision arithmetics) the 
remedy is easy. One may take, say, the A^ = 3 toy metric of Ref. [2], 



ef)(a) = 



2 2ga 







2 ga a + 1 g 

a+2 







9 



2a+l 



and represent the triplet of eigenvalues Pj{g), j = 1, 2, 3 (as sampled in Figure 
Nr. 3 of loc. cit. at a = 1) in logarithmic scale yielding, say, the adapted 
= 1 secular equation 



a 



det 



2 - e-"'^-20 
2g 2 




^g 

9 





9 



i.e., the non-polynomial version of our eigenvalue problem. 
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In a test run the numerical analysis of this equation reproduced the results 
given in Table Nr. 1 of loc. cit.. 

On this basis one may expect that the key problem brought by the rescal- 
ing appears tractable by the MAPLE-based numerical software. The main 
gain came with the substantial extension of the feasibility of the graphical 
determinations of the parameter-dependence of the physical domains "Pg^"* (a) 
with the growth of N. The characteristic illustration is offered by Figures [1] 
and [2] which clearly demonstrate the emergence of an obvious pattern which 
was not accessible without rescaling. 




Figure 1: Seven eigenvalues p = p{g) of metric Q^J\l 
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Figure 2: Nine eigenvalues p = p{g) of metric 'df\l 



5 Summary 



From the point of view of Quantum Mechanics it is rather unfortunate that 
for a given Hamiltonian H the specification of the metric Q{H) prescribed by 
Dieudonne Eq. ([1]) is ambiguous [9]. In this context, our series of papers [1] - 
[7] has been devoted to the constructive study of the one-to-many mappings 
H — )■ Q[H). In essence, we offered there a new methodical recipe of a 
systematic suppression of the ambiguity of the menu of eligible Q{H)s. 

In our present continuation and extension of this paper we decided to 
explain the symbolic-manipulation aspects of such a recipe in more detail. 
Emphasizing that such a problem would be hardly tractable and/or solvable 
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without an essential interaction between the abstract quantum theory and 
the symbohc-manipulation techniques and algebraic constructionis assisted 
by contemporary computers. 

One of byproducts of such an interaction between methods has been 
shown to lie in the possible amendments of the numerical aspects of the 
necessary simultaneous analysis and study of physical domains and of the 
other properties of both of the physics-representing operators H and B. On 
the basis of these results one can conclude that the confirmation of feasibility 
of a methodical symbiosis between algebra and analysis seems to be able to 
contribute to the contemporary quick growth of popularity of the practical 
phenomenological applications of quantum models with nontrivial metrics, 
e.g., in optics (l8j . 



Appendix A. A compact review of Quantum 
Mechanics using nontrivial metrics Q ^ I 

In the spirit of certain textbooks on Quantum Mechanics [19] the Dirac's ket 
vectors may be interpreted as describing different quantum systems 

the nature of which may vary with the other characteristics of the system 
in question. Thus, one can start from the most common single-particle (or, 
in some applications, single quasi-particle) quantum models where, roughly 
speaking, there exists an operator Q of the particle position with eigenvalues 
g G M and eigenkets One then introduces the standard (often called 
Dirac's) Hermitian conjugation (i.e., transposition plus complex conjugation) 

^{Dirac) . ^ _ (11) 

In this manner one represents the kets of states by the concrete square- 
integrable functions ip{q) = (qlip) (of the coordinate q) called wave functions. 
In this context the discrete, finite-dimensional models exemplified by Eq. ([3]) 
may be perceived as the most convenient testing ground of the methods. 

In all of the similar classes of examples the vector space of states V be- 
comes endowed with the most common inner product 

/C30 
ra{q)Mq) dq (12) 
-oo 

with, possibly, the integration replaced by the infinite or finite summation. 
Thus, we may (and usually do) set V = i^^(IR), etc. 

Without any real danger of misunderstanding we may speak here about 
the "friendly" Hilbert space of states H^^^ = V, calling the variable q in 
Eq. f|T2|) "the coordinate" . In parallel we usually perform a maximally con- 
venient choice of the Hamiltonian H based on the so called principle of corre- 
spondence which "dictates" us to split H = T + V where the general interac- 
tion operator V is represented, say, by a kernel V{q, q') when acting upon the 
wave functions. In practice, this kernel is most often chosen as proportional 
to the Dirac's delta-function so that V becomes an elementary multiplicative 
operator V = Viocai = ^{l)- Similarly, the most popular and preferred form 
of the "kinetic energy" T is a differential operator, say, T = Ti^cai = —d? jdq^ 
in single dimension and suitable units. 
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A word of warning emerges when we perform a Fourier transformation 
in H^^^ so that the variable q becomes replaced by p (= momentum). One 
should rather denote the latter, "Fourier-image" space by the slightly differ- 
ent symbol T-L^^\ therefore (with the superscript still abbreviating "physical" 
|11|). Paradoxically, after this change of frame the kinetic operator Tiocai be- 
comes multiplicative while Viocai becomes strongly non-local in momenta. 

The truly deep change of the paradigm comes with the models where the 
necessity of the observability of the coordinate q is abandoned completely. 
One may still start from the vector space of kets V and endow it with the 
Hilbert-space structure via the similarly looking inner product 

m^b) = j^ra{s)i^b{s)ds. (13) 

The non-existence of the position operator Q changes the physics of course. 
The key point is that we lose the one-to-one correspondence between the 
integration path C and the spectrum M of any coordinate- mimicking operator. 
The physics-independent optional variable s remains purely formal. 

In such a setting our choice of the physical observables remains unre- 
stricted and is entirely arbitrary. It need not be related to any usual classical 
system, either. For illustration one might recall the pedagogically motivated 
paper PU]. In a slightly provocative demonstration of the abstract nature 
of quantum theory the variable s in Eq. (fT3|) has been chosen there as an 
observable "time" of a hypothetical "quantum clock" system. 

Definition ( !T3|) of the inner product in 'H*^'^) or in is the theoretical 
framework within which the traditional quantum mechanics works, reducing 
the full power of the theory to something acceptable via analogies with clas- 
sical physics. Still, one need not move too far. Even in Ref. [T], for example, 
we did not treat the variable s in Eq. f|T3|) as a purely formal quantity, having 
preserved at least a part of its relation to the position q. In particular, we 
stayed in the middle of the path towards abstraction and worked still with 
the usual one- dimensional ordinary differential Schrodinger equation 

--ni^n{s) + V{s)iJn{s)=Eni^n{s), n = 0, 1, . . . (14) 

where the key nonstandard features can be seen in 

• the admissibility of the complex potentials sampled by the power-law- 
anharmonic family V{s) = — (is)^"'"'' of Ref. [21] giving the "standard" 
real and discrete bound-state spectra at any S > (cf. the proofs in 

my, 

• the admissibility of the replacement of the usual real line of s by a 
suitable complex curve C = C(M) which may even be, in principle, 
living on a complicated multisheeted Riemann surface [23] : 

• in the possibility of a systematic study of its discrete analogues and 
simplifications. 

Naturally, one partially leaves the more or less safe guidance offered by the 
principle of correspondence. Just a partial revitalization of this guidance is 
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possible (cf., e.g., a nice example-based discussion of this point in Ref. [21]). 
A partial reward for this loss can be sought in a new tractability of some 
traditionally contradictory quantization problems (say, when working, say, 
with the relativistic Klein-Gordon equations [23]). 

The main theoretical difficulty consists in the vast ambiguity of the nec- 
essary appropriate generalization of Eq. (fTTl) . The general recipe (explained 
already in [Qj or, more explicitly, in [11]) is metric-dependent and reads 

^^.ene^O^I^^ ^ U\ := . (15) 

This means that using the language of wave functions ipi^) with s G C we 
must replace the most common single-integral definition (fT3|) of the inner 
product in the original "friendly" Hilbert space Ti^^^ by the more sophisti- 
cated double-integral formula 

{{^a\iJb)= [ [ ra{s)&{s,s')Ms)dsds'. (16) 

Jc Jc 

In terms of an integral-operator-kernel representation 0(s, s') of our abstract 
metric operator = > this recipe defines a more sophisticated inner 
product which converts the ket-vector space V into the metric-dependent 
(and physics- representing, "standard" [TT]) Hilbert space Ti^^^ . 
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